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ABSTRACT 

Three-dimensional hydrodynamic simulations exploring the first 18 hours of the 2010 January 28 
outburst of the recurrent nova U Scorpii have been performed. Special emphasis was placed on 
capturing the enormous range in spatial scales in the blast. The pre-explosion system conditions 
included the secondary star and a flared accretion disk. These conditions can have a profound influence 
on the evolving blast wave. The blast itself is shadowed by the secondary star, which itself gives rise to 
a low-temperature bow-shock. The accretion disk is completely destroyed in the explosion. A model 
with a disk gas density of 10^^ cm^"^ produced a blast wave that is coUimated and with clear bipolar 
structures, including a bipolar X-ray emitting shell. The degree of coUimation depends on the initial 
mass of ejecta, energy of explosion, and circumstellar gas density distribution. It is most pronounced 
for a model with the lowest explosion energy (10'*'^ erg) and mass of ejecta (IO^^A/q). The X-ray 
luminosities of three of six models computed are close to, but consistent with, an upper limit to the 
early blast X-ray emission obtained by the Swift satellite, the X-ray luminosity being larger for higher 
circumstellar gas density and higher ejecta mass. The latter consideration, together with estimates of 
the blast energy from previous outbursts, suggests that the mass of ejecta in the 2010 outburst was 
not larger than 10~^ Mq. 

Subject headings: stars: novae — stars: individual (U Sco) — X-rays: binaries — shock waves — 
methods: numerical 



1. INTRODUCTION 

U Scorpii was discovered to have entered its latest 
outburst on 2010 January 28 by amateur astronomer 
B.G. Harris ( Schaefer et al .112 01 Of) . It is one of the best 
observed of th e 10 recurrent novae known to date (e.g. 
ISchaefeiil2010l) . with a fairly constant recurrent cycle of 
8-12 yr and a previous outburst in 1999. It comprises a 
white dwarf of mass 1.55 ± O.24M0 and late- type star of 
m ass 0.88 ± 0.17Mp^ in an orbit with a period of 1.23d 

f)hnston fc Kulkarnil 119921: iThoroughgood et~all 120011: 
laefer 1990: iSchaefer fc Ringwaldl I1995D . iSchaeferi 
IIQ) adopts a median spectral type of the secondary 
from estimates culled from the literature of G5 IV. 

The outburst mechanism for recurrent novae is ther- 
monuclear runaway on the white dwarf triggered by the 
accre t ed ma ss exceeding a critical limit (jStarrfield et al.l 
119851 I1988D . The relatively short 10 year cycle of 
U Sco and the fast evolution of its outbursts im- 
plies the white dwarf is c lose t o the Chandrasekha r 
limit (e.g. IStarrfield etldl 119881: iHachisu et all |2000D . 
Ejecta from previous outbursts were also found to 
be He-rich (|W illiams^e LaL| [ 19811; iBarlow et all 119811: 
lAnupama fc De wangan l2000t though this was con- 
tested bv lliiimal 120021 ). further fueling discussion of 
U Sc o as a possible progenitor of a Typ e la s upernova 
(e^ lStarrfield et al.lfl988L IHachisu et al]ll999l : ISchaefed 



of prompt X-ray observations. RS Oph was a particu- 
larly bright X-ray source: the explosion occurs within 
the dense wind of the secondary red giant star, result- 
ing in a ~ 10^-10^ K blast wave whose evolutio n was 
followed by RXTE, Sw i ft and Chandra ([S okolosk i et al. 



2006; 'Bod e et al.l |2006; 'Nel son et afl [20081 : IDrake et al 
2009 , 1. Detailed SD hydr odynamic modelling of the blast 



bv [Orlando et al.l (|2009[ ) provided estimates of the ex- 
plosion energy, ejecta mass, and confirmed the presence 
of enhanced gas density in the equatorial plane of the 
system. 

Spectroscopy of the 2010 outburst of U Sco found 
Ha and Up H Balmer line profiles broadened by 
7600 km s~^ full- width at half- maximum (|Anupamal 
120101) ■ and bv 11000 km s^^ full- width at zero inten- 
sity ([Aral et al.ll20100 , consistent with those observed in 
earlier outbursts Ce.g. IBarlow et al 119811: [Munari et al.l 



20101 see, howev er, the evolutionary co i isider ations of 



Sarnaet al.ll2006D : IThoroughgood et al.l (|200lD charac 



terised it as "the best Type la supernova progenitor cur- 
rently known" . 

The 2006 explosion of the recurrent nova RS Oph pro- 
vided a cogent demonstration of the diagnostic power 



119991: lAnupama fc Dewanganll2000l : IIiiimall2002D . If par- 
tially converted to thermal energy, analogous to the ex- 
plosion of RS Oph into the dense wind of the secondary 
star, such kinetic energy would result in heating to X-ray 
emitting temperatures. While the secondary of U Sco 
lacks the cool, dense wind of the RS Oph secondary that 
provided the medium for copious X-ray production, the 
possibility of the blast generating significant X-rays from 
shocked ejecta or shocked gas in its accretion disk or am- 
bient circumbinary material remained. 

The early outburst was observed by the Swift, 
RXT E and Integral sate l lites within a d ay of discov- 
ery (|Schlegel et al.l l2010t lManousakis"er al. 201 Q) but 
all fai led to detect any X-ray emission. ISchlegel et al.l 
(|2010t) placed the most stringent upper limit of Lx < 
1.1 X 10"" erg s-icm-2 (90% confidence) in the 0.3-10 



keV band based on data obtained on 2010 Jan 29-30. 

Here we describe detailed hydrodynamic simulations 
of the U Sc o explosion sim i lar to those undertaken for 
RS Oph by lOrlando et all (|2009D . We investigate the 
effects of the accretion disk and close secondary compan- 
ion on the explosion, and examine possible constraints 
the lack of observed X-rays might provide. 

2. HYDRODYNAMIC MODELLING 

The 3-D hydr odynamic mod e l adop ted here is sim- 
ilar to that of lOrlando et~all ()2009l) . The calcula- 
tions were performed using FLASH, an adaptive mesh 
refine ment multiphysics code for astrophysical plas- 
mas (jFryxell et al.l |2000[ ). The hydrodynamic equa- 
tions for compressible gas dynamics are solved using 
the FLASH im plementation of the piecew ice-parabolic 
method ( PPM: IColella fc Woodwaldl [1981 . Unlike the 



models of lOrlando et al.l (|2009D computed for the case of 
RS Oph, for the very early phase of the blast studied here 
we do not take into account radiative losses and thermal 
conduction, assuming that the remnant is in its adiabatic 
expansion phase for the whole evolution considered. This 
assumption is valid if the time covered by our simulations 
iend is smaller than the transition time from adiabatic to 
radiative phase for an expanding spherical blast, defined 
as (e.g. jBlond in et al. 1998; Pctruk 2005), 



itr = 0.522 E^l^'^ 



-9/17 



(1) 



where E = E^va the energy of the explosion and n^^A is 
the particle number density of the ambient medium. As 
discussed below, it turns out that icnd ■^ ^tr in our set of 
simulations and our assumption is valid. 

For the system paramet ers, we adopt the values of 
iThoroughgood et al.l (|200lD : these are listed in Table [H 
together with the parameters adopted for the blast mod- 
els investigated here. We included a diffuse, puffed- 
up disk-like distribution around the system which is 
analogous to the gas density enhancement found in 
the numerical model of the RS Oph sys t em im me- 
diat ely prior to outburst by IWalder et al.l (|2008l see 
also iMastrodemos &: Morrij I1999D and to the equato- 
rial gas distribution adopted by lOrlando et al.l ()2009l ). 
We also included a flared, uniform density gas accretion 
disk around the wh ite dwarf, adopting the shape from 
iHachisu et al.l (J2000D . We included the secondary star as 
an impenetrable body. A close-up of this initial pre-blast 
configuration is illustrated in Figure [1] 

As an initial c ondition, we assumed a spherical Sedov- 
type blast wave ()Sedovlll959D originating from the ther- 
monuclear explosion centered on the white dwarf, with 
radius rg = 4 x 10^ km and total energy Eq. This 
energy was partitioned such that 1/4 was contained in 
thermal e nergy, and 3/4 in kinetic energy. Solar abun- 
dances of iGrevesse fc Sauval (|1998D (GS) were assumed 
for the disk, wind and equatorial ambient gas, while the 
ejecta metallicity was assumed enhanced by a factor of 
ten. This latter choice was guided by observations of 
He-rich ejecta in previous o utbur st, as noted in Sect. [U 
as well as the iDrake et al.l (|2009D high-resolution X-ray 
spectroscopic study of the 2006 RS Oph blast that found 
evidence for metal-rich ejecta. Since radiative losses were 
not important for the early blast evolution, the choice of 
abundances is only relevant for computation of the emit- 
ted X-ray intensity of the blast. 




Fig. 1. — A rendition of the U Sco binary system model initial 
conditions. The white sphere represents the pre-outburst, initial 
radius of the Sedov-type blast on the white dwarf. 



Six models were computed to explore the effects of a 
different accretion disk and circumstellar gas density, ex- 
plosion energy and total ejecta mass on the outcome of 
the blast. Direct estimates of the ejecta n iass from obser 
vations of the 1979 and 1999 outbursts fWilliams et al 
1981'i Anupa ma fc D ewangan 2000; Evans ct al. 200lf 
find Me.j ^ 1O~^M0. Expansion velocities of 






5000 km s~^ from optical spectroscopy correspond to 
a kinetic energy of 2.5 x 10*^ erg, which can be taken 
as a lower limit to the explosion energy. This is consis- 
tent with estimates of the integrat ed optical outpu t from 
earlier outbursts of 10'''^ erg fe.g. IWebbink et all [1987). 
We therefore adopted fiducial values of Eq = lO^** erg 
and Afej — 1O~^M0 and explored the effects of lower 
and higher values of each by an order of magnitude. 
For the accretion disk we assumed a uniform density 
of Tid = 10^^ cm~^ and investigated a further model 
with Ud = 10^^ cm~'^. The wind and equatorial den- 
sity enhancement are included for completeness. While 
this circumstellar gas is not expected to have a profound 
infiuence on the dynamics and evolution of the blast, it 
does influence the shock structure and the predicted X- 
ray luminosity. We probed the effects of increasing and 
decreasing the circumstellar gas density by a factor of 10. 

The hydrodynamic equations were solved in one quad- 
rant of the 3-D spatial domain and the coordinate system 
was oriented in such a way that both the white dwarf and 
the donor star lie on the x axis. The donor star was at the 
origin of the coordinate system, (a;o,yoj2o) = (0,0,0), 
and the white dwarf was located at a; = 0.03 AU, i.e. 
the system orbital separation. The extents of spatial 
domains employed for the different models are listed in 
Table [1] together with other relevant parameters, and de- 
pended on the circumstellar gas density; two models that 
investigated lower densities for which expansion rates 
were greater required larger domains. 

The explosion was followed for 18 hours in order to 
explore the early phase X-ray emission and capture the 
details of the effects of the explosion environment on the 
blast evolution. A major challenge in modelling the ex- 
plosion of U Sco is the very small scale of the system com- 
pared with the size of the rapidly expanding blast wave. 
The separation of the two stars is only 0.03 AU — a dis- 
tance traversed by a blast wave moving at 5000 km s~^ in 



TABLE 1 

Adopted parameters and initial conditions for the hydrodynamic models of the 2010 U Sco explosion 



Parameter 








Value 






Secondary star radius 








Rb = 2.1R0 






Binary separation 








a = 6.5ij0 






Inclination 








i = 83 deg 






Model abbreviation 




E44M-7D13 


E44M-7D15 


E43M-8D15 E44M-7D15HW 


E44M-7D15LW E45M-6D15LW 


Explosion energy Eq (erg) 




10" 


10« 


IQAi IQiA 


10« 


10^!' 


Ejecta mass M^j (Mq) 




10-7 


10-7 


10-8 10-7 


10-7 


10-'' 


Accn. disk density n^ (cm' 


-■') 


10" 


1015 


XQis iQis 


IQis 


1015 


Sec. wind density n„ (cm~ 


--') 


10« 


106 


10'5 107 


iqs 


10^ 


Eq. enhancement rieq (cm" 


-■') 


10* 


lOS 


108 109 


107 


107 


X-ray luminosity Lx (erg i 


5-^) 


1 X 10^3 


5 X 10^2 


3 X 10^1 4 X 10^^ 


5 X 10^1 


7 X 1033 


Swift XRTLx (ergs-i) 








< 2 X 10^3 (90% confidence) 






Spatial domain 






-6.6 < 


X < 6.6 AU 


-13.2 < 


X < 13.2 AU 








0<y 


< 6.6 AU 


0<y 


< 13.2 AU 








0< z 


< 6.6 AU 


0<z 


< 13.2 AU 


AMR max. resolution 








4 X 108 cm (2.7 x 10-^ AU) 






Time covered 








0-18 hours 







less than 20 minutes — and the secondary star subtends 
a sufRciently large solid angle that it has a significant 
effect on the blast wave evolution. To capture this range 
of scales, the models explored here employed 16 nested 
levels of adaptive mesh refinement, with resolution in- 
creasing twice at each refinement level. This grid config- 
uration yields an effective resolution of « 4 x 10'' km at 
the finest level, corresponding to « 100 grid points per 
initial radius of the blast. 

Given the parameters in Table [TJ for five of our six 
models ttr > 3 days from Eq. [H and for the sixth (E44M- 
7D15HW) ttr — 2.3 days. These times are much larger 
than the time covered by our simulations. Our modeled 
remnants therefore never enter the radiative phase. 

3. RESULTS AND DISCUSSION 
3.1. Temperature and density structure 

The U Sco models all evolved quite different ly to those 
presented for RS Oph bv lOrlando et all ()2009[ ). As noted 
in §1, this difference is expected because U Sco lacks 
the dense wind of RS Oph for propagation of an X-ray 
emitting blast wave. 

The 18 hour post-blast gas density and temperature 
distributions in the x-z plane bisecting the system (the 
plane of the orbital axis; the accretion disk is edge on) 
are illustrated in Figure [2] All the models exhibit con- 
spicuous departure from spherical symmetry. The main 
feature of the gas density distribution is a large cavity 
to the left of the figures (negative x direction) which is 
a blast wave "shadow" cast by the secondary star. This 
result reinforces the need to perform simulations with 
sufficient spatial resolution to capture the details of the 
pre-blast system configuration. 

In all models the accretion disk is completely destroyed 
by the blast. This is not surprising considering the grav- 
itational binding energy of the disk in our models is 3- 
4 orders of magnitude lower than the explosion energy: 
disk destruction is assured unless the di sk gas density 
is sub stantiall y higher than 10^^ cm"^. IWorters et al.l 
(|201Ci see also lMunari et al.ll2010( ) report a renewed op- 
tical flickering with amplitude "0.2 mag over the course 
of an hour" on 2010 February 5. They interpret this 
re-establishment of accretion, implying a limit of 7 days 



on the build-up of an angular momentum shedding ac- 
cretion disk. In our simulations, the disk volume is 
Vdisk « 2 X 10"^^ cm'' and its mass ranges between 
2.2 X 10-i°Mq (model E44M-7D13) and 2.2 x lO-^'Mg 
(all other models). The build up of the disk in 7 days 
would imply a rate of mass loss from the companion star 
ranging betwee n 10~^ and 10~^Mp) y r-^, consistent with 
the estimate of iHachisu et al.l (|2000( ). We note, however, 
that any renewed disk build-up would need to overcome 
the dynamic pressure of the radiatively-driven outflow 
from the evolving central object. 

Model E44M-7D13, with lower accretion disk gas den- 
sity, shows only a very small effect of the disk in a slight 
inhibition of the blast wave progress in the positive x di- 
rection. In contrast, the higher disk gas density in all the 
other models has quite a profound effect on the blast evo- 
lution. Both the blast and ejecta are strongly collimated 
in polar directions, the collimation being more prominent 
for lower explosion energy (E43M-8D15) as might be ex- 
pected. Figure [3] shows, as an example, the collimation 
of ejec ta 18 hours after the out burst for the model E44M- 
7D15. iKato fc Hachisul ()2003l ) interpreted "horned" op- 
tical line profiles observ ed in previous U Sco outburst 
(e.g. iMunari et al.l 119990 as evidence for "jet-shaped" 
outflows. Similar blast collimation was both predicted 
by hydrodynamic simulations of the early R S Oph blast 
by Waldcr et al. (2008) and Orlan do et all dfOOQ), and 
deduced through X -ray spectroscopy obtained 13.6 days 
after the outburst (jDrake et al.l 120091 ). The simulations 
indicated the collimation was caused by an equatorial 
disk-like gas distribution around the white dwarf. High 
spatial resolution observations of the blast at various 
times also found o ptical, infrared, radio and X-ray colli- 
mation signatures (lO^Brien et al.ll2006l: iBode et al.ll2007l: 
IChesneau et al.ll2007t iLuna et al.b'ooa) . 

Evidence for a more strong ly collimated synchrotron 
radio jet was a lso fo und by iRupen et al.l (J2008) and 
ISokoloski et al.l (|2008D . The latter authors concluded 
that such structure is too highly collimated to be ex- 
plained by the influence of circumbinary material on 
the explosion. The calculations presented here and by 
I Orlando et al.l (I2009D bring such conclusions into ques- 
tion: the nature of the explosion environment is cru- 
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Fig. 2. — Colour-coded cross-section images of the gas density distribution (top) and temperature (bottom) 18 hours after the U Sco 
blast, corresponding to the end of the hydrodynamic simulations. The secondary star is at the origin, and the white dwarf to the right. 
The different panels illustrate models with accretion disk density n^ = 10^'^ cm~'^ and an explosion energy of Eq = 10^^ erg (E44M-7D13; 
Left); nd = lO^'^ cm-^ and Eo = 10** erg (E44M-7D15; left-center); rij = lO^'^ cm-^ and Eo = 10*^ erg (E43M-8D15; right-center). Far 
right panels show a similar model to the left-center panels, except with CSM gas densities increased by a factor of ten (E44M-7D15HW). 
Inset panels show the blast structure closer to the system origin. The secondary star is shown and is responsible for the strong "shadowing" 
of the blast wave. The accretion disk is completely destroyed by the blast in all simulations. The white dotted contour encloses the ejecta. 
The red solid contour denotes the regions with plasma temperature T > 1 X 10'' K. 

enhancements in the ejecta. 

The outer regions of the blast are characterised by 
a layer of shocked circumstellar medium (CSM) with 
temperature « 10^ K and a thin expanding surface of 
shocked ejecta with temperature in the range [10^ — 
10*] K (see Figure ^. The CSM gives rise to a re- 
flected shock that heats the ejecta, and is strongest for 
model E44M-7D15HW in which the CSM density is high- 
est. The secondary star also gives rise to a promi- 
nent "bow shock", though in this structure the gas is 
shocked to a temperature of only a few 10^ K and does 
not give rise to any conspicuous observable X-ray emis- 
sion. It is worth noting that Fig. [5] reports the pro- 
ton temperature set by coUisionless heating at the shock 
front. Given the short timescale considered here, pro- 
tons and electrons are likely not thermalized at the for- 
ward shock, so that the ratio of the electron to pro- 
ton temperat ure at the blast wave is expected to be 
(Te/Tp) < 1. iGhavamian et al.l ()2007D derived a phys- 
ical model for the heating of electrons and protons in 
strong shocks that predicts a relationship (Te/Tp) ex v^ , 
where vs is the shock velocity. Considering the velocity 
of the forward shock in our s imulations is of the orde r of 
10000 km s'\ from Fig. 2 of :Ghavamian et all (|2007D we 
derive (To/Tp) < 0.1. Thus, we expect that the electron 
temperature of the shocked CSM is « 10® K. X-rays are 
instead naturally dominated by the higher density parts 
of the shocked ejecta which are expected to equilibrate 
quickly. In the following, we therefore assumed electron- 
proton temperature equilibration in shocked ejecta. 




Fig. 3. — Three-dimensional rendering of the ejecta distribution 
18 hours after the blast for the model E44M-7D15. The plane of 
the orbit of the central binary system lies on the (x, y) plane. 

cial for the morphology of the evolving blast and ejecta. 
While the models presented here do not produce obvious 
narrow jets of material, we note that radio synchrotron 
emission also depends strongly on the magnetic field 
stren gth, gradient and orientation (e.g., I Orlando et alj 
l2007f ) and need not be co-spatial with obvious density 




Fig. 4. — X-ray images (normalized to the maximum of each panel) projected along the line of sight of the blast after 18 hours of 
evolution, corresponding to the density and temperature distributions illustrated in Figure[2] Upper panels are rendered with a logarithmic 
surface brightness scale and lower panels with a linear scale. The plane of the orbit is inclined by 7 degrees to the line-of-sight. Both the 
secondary subgiant and the white dwarf lie on the x axis. Most of the X-ray emission originates from shocked ejecta. 



3.2. X-ray luminosity and morphology 

X-ray emission from the blast was synthcsised from 
the model res u lts us ing the methodology described by 
I Orlando et al.l (|2009f) . The synthesis includes thermal 
broadening of emission lines, the Doppler shift of lines 
due to the component of plasma velocity along the line- 
of-sight, and photoelectric absorption by the interstel- 
lar medium (ISM), CSM, and ejecta. The absorption 
by the ISM is calculated as suming a column de nsity 
A^H = 1.4 X 10^1 cm-2 (e.g.. iKahabka et allfTool and 
is consistent w ith a distance to U Sco of 12 kpc, e.g. 
ISchaefeiil201Cl( l: the local absorption is calculated as due 
to shocked CSM (assuming GS abundances) and ejecta 
(GS abundances xlO). The X-ray emission from the 
blast at 18 hours is illustrated for four representative 
models in Figure U) The absorption by the ejecta con- 
tributes to a small difference between the predicted X- 
ray luminosities of the different models. Most of the 
emission arises from the high-temperature shocked ejecta 
forming the outer "shell" of the blast wave where we as- 
sumed electron-proton temperature equilibration. Model 
E44M-7D13 exhibits significant equatorial brightening 
that arises from a ring of higher density ejecta that was 
partially impeded by the accretion disk. This structure 
is not present in models E44M-7D15, E43M-8D15 and 
E44M-7D15IIW, with a higher density disk. Instead, 
these latter models are dominated by bipolar cusps that 
appear either limb-brightcncd (E44M-7D15), or forming 
ring-like structures (E43M-8D15 and E44M-7D15HW), 
in the projected image, accompanied by limited equato- 
rial emission. 

The X-ray luminosities of the simulations in the [0.6- 
12.4] keV band are listed in Table [T] and vary between 
Lx = 3 X 10^1 erg s'^ (E43M-7D15) and 7 x 10^3 erg s'^ 



(E45M-6D15LW). The X-ray luminosity is larger for 
higher ejecta mass, explosion energy, and circumstellar 
gas density. The X-ray flux upper limit of Fx < 1.1 x 
10"^'^ erg s~^c m~^ obtained at a s imilar epoch to our 
simulations by ISchlegel et al.l ()2010t ) corresponds to an 
X-ray luminosity of Lx = 2 x 10'^^ erg s~^ for a distance 
of 12 kpc. The fact that the models with M^j = 10"'^Mq 
predict X-ray luminosities close to the Swift upper limit 
(except run E44M-7D15LW whose luminosity is much 
lower), while the model with Afcj = lO^^M© and lower 
circumstellar density (E45M-6D15LW) exceeds it, sug- 
gests that the mass of ejecta in the 2010 outburst was 
not significantly larger than 1O~^M0. This is consis- 
tent with esti mates of previous outbursts in 1979 and 
1^99j[Williams,.etalJ[1981|; lAnupama fc Dewanganll2000l: 
I Evans et al .1120011 ). While we cannot rule out with cer- 
tainty the circumstellar density being even lower than 
assumed in the latter model with Afcj — 1O~^M0, which 
could reduce the predicted X-ray luminosity and recon- 
cile it with the Swift upper limit, such a large mass re- 
quires the larger assumed explosion energy of 10'*^ erg to 
take the ejecta out of the gravitational well of the system. 
This exceeds by an order of magnitude the total explo- 
sion energ^f_ofJ^O^ergestimated for previous outbursts 
by (IWebbink et allflOSl . 

The model with a higher wind and equatorial gas den- 
sity (E44M-7D15IIW) has a stronger X-ray contribution 
from ejecta heated in the reflected shock that exceeds 
the observed X-ray upper limit. This is consistent with 
our assumption of a rather low circumstellar gas density, 
commensurate with the evolutionary state of the U Sco 
component stars, and does not alter the conclusion that 
the ejected mass was probably not larger than 1O^^M0. 
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